Kinetics of the disordered Bose gas with collisions 
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We present a semi-phenomenological approach to the kinetics of the disordered interacting Bose 
gas, based on the Boltzmann transport equation. The theory may serve as a unifying framework 
for studying questions of dynamics of the expanding Bose gas at different stages of the expansion. 
We show that the transport theory allows to straightforwardly reproduce and generalize a number 
of results previously obtained from microscopic models in different formalisms. Based on estimates 
for the inter-particle scattering rates, we discuss the relevance of interaction effects for the localiza- 
tion problem in the interacting disordered Bose gas. We argue that, if the number of particles is 
large enough, the size of the expanding cloud may exceed the localization length. We describe the 
spreading of the wave packet in this regime as collision-induced diffusion and compare the obtained 
rate of expansion to known results on subdiffusive spreading in non-linear disordered lattices. 

PACS numbers: 71. 10. Ay, 71.10.Pm, 75.40.Cx 



m 
o 

(N 



<Z3 

wo. 



i 

C 

O 
o 



> 

00 
(N 
O 

o 

O 

m 



X 



Introduction: Studies of localization effects in disor- 
dered Bose gases are at the forefront of research on 
cold atomic gases [ll-Hoj. In a typical setup, the gas is 
released from a trap and expands under the influence 
of a disorder potential. The evolution of the expand- 
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ing gas cloud is often divided into two stages 
the first stage, the density is comparatively high, the 
interparticle-interaction dominates, and disorder-related 
effects are neglected. At the following stage the density of 
the cloud drops so that the average kinetic energy (per 
particle) becomes much larger than the interaction en- 
ergy. For this second stage, the phenomenon of Ander- 
son localization for non-interacting particles is studied. 
As we discuss in this paper, this picture is not complete 
due to the effects induced by interparticle-collisions. 

In Ref . Il2| we raised the question how the presence of 
interactions influences the expansion of a matter wave- 
packet in a disorder potential. During the course of the 
expansion, the energy stored in an interaction-induced 
Hartree-potential is gradually converted into kinetic en- 
ergy, and the combined effect of disorder and interac- 
tions leads to interesting dynamics. This study as well 
as the related work of Ref. [l3| were based on the Gross- 
Pitaevskii equation (GPE), thereby assuming a large oc- 
cupation of participating modes that allows to describe 
the Bose gas as a non-linear medium. Initially, only the 
collision-less limit was considered, while the effect of col- 



lisions has been studied only recently [14 1 



In this work, building upon insights gained from the 
microscopic analysis of the GPE, the appropriate ki- 
netic equation for the disordered interacting Bose gas 
is obtained from a semi-phenomenological Boltzmann- 
equation approach. The underlying physics is essentially 
classical: Throughout the paper we assume that the in- 
teraction energy is sufficiently small compared to the ki- 



netic energy so that one may neglect off-diagonal com- 
ponents of the Bogoliubov Hamiltonian. 

The classical theory of nonlinear diffusion based on 
the Boltzmann equation can serve as a starting point 
for a qualitative analysis of the effect of interparticle- 
interactions on localization. There are three main effects. 
First, the time-dependent interaction- induced potential 
and the inter-particle collisions can both lead to dephas- 
ing. In particular, the inter-particle collision rate is a 
lower bound for the dephasing rate when considering ef- 
fects of localization. Secondly, for an out-of equilibrium 
distribution, collisions lead to a particle-flux down in en- 
ergies. Since particles with a low energy are much more 
likely to localize, this process may serve as a seed for 
localization. Thirdly, the wave packet may continue to 
spread as a result of interparticle collisions despite the 
fact that single particle states are localized. The find- 
ings can be important for experiments on disordered Bose 
gases as soon as interactions are not completely negligible 
on the stage of localization. 

In the following we derive an effective nonlinear dif- 
fusion equation starting from the semi-phenomenological 
Boltzmann equation. We consider systems for which dis- 
order scattering is the most effective scattering mecha- 
nism, while inter-particle collisions are sufficiently infre- 
quent. In the resulting equation, the interaction-induced 
collision integral simplifies considerably since in the dif- 
fusive limit the distribution functions may be taken as 
isotropic. This allows us to provide analytic expres- 
sions for the kernel of the collision integral in the two 
and three-dimensional cases. We point out that for a 
two-dimensional gas the mean squared radius is a lin- 
ear function of time, whenever the impurity scattering 
can be modeled by a white noise random potential. This 
result can be used as a reference point for numerical sim- 
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ulations. Next, we make contact both with our recent 
work based on the GPE in the limit of large occupation 
numbers 1J| and with the work on thermalization of dis- 
ordered Bose gas for low occupation numbers, Ref. [15. 
Finally, we discuss the spreading of the wave-packet re- 
sulting from interparticle-collisions. We predict that if 
the number of particles is large enough, the size of the 
expanding cloud may exceed the localization length. 

Boltzmann equation: We consider the Boltzmann 
equation in the form 



d t f p (v,t) +vV/ p (r,i) - V£/(r)V p /p(r,i) 

= C*22 [/] + Ci m p [/] • 



(1) 



Here, 



p/m is the velocity, the distribution func- 



tion that is normalized as n = J(dp)f p , where 
the density of particles and we use the abbreviation 
(dp) = dp/(27r) d in d dimensions. The potential U in 
general comprises two parts, a smooth classical exter- 
nal potential, and a self-consistent Hartree-like constri- 
bution, #(r,i) = 2g ( in(r 1 t) [H, [l]|. The coupling con- 
stant in three dimensions equals 53 = 47ra s /m, where a s 
is the scattering length; in two dimensions one should set 
.92 = .93 /<iz, where a z is the length which characterizes 
the transverse confinement. The interpretation of the 
above equation is a statistical one: We consider quanti- 
ties that are averaged over many realizations of disorder. 
This is why the density n, for example, is a smooth quan- 
tity. The equation should be supplemented with initial 
conditions, alternatively one may add a source term to 
the right hand side of the equation. 

Depending on the experimental setup, specifics of the 
system like the finite correlation length of the disorder in 
a speckle potential need to be taken into account lil 17|. 
For example, the smooth external potential U may incor- 
porate the Fourier components of the speckle potential 
that are much smaller than A -1 (here A is the typical 
wavelength of the Bose-gas particles), while the shorter 
components are comprised into the impurity collision in- 
tegral Cimp on the right-hand side of the above equation. 
For the sake of simplicity, we will consider a simple re- 
laxation type approximation, for which 



1 



(/ P (r) - </p(r)>, 



(2) 



and (■•■)„ denotes an averaging over the solid angle of 
the momentum. r p is the transport scattering time. The 
collision integral C22 for the interacting Bose gas with a 
short-range pseudo-potential reads (Hf , 

C 22[f] = (2^)2rf / d P2 d P3^P4 S(p + p 2 - p 3 - p 4 ) 



xS(e p + e 



PI °P3 C P4 



^WL (3) 



where 



HI] = (1 + /p + /pJ/pa/p* " /p/ P2 (1 + /ps + / P J (4) 



is a combination of distribution functions characteristic 
for the in- and out-scattering processes. For the sake 
of brevity, we suppressed the space arguments for the 
distribution functions. It is worth mentioning that in the 
limit of large occupation numbers, i.e., when the GPE as 
a classical field equation is valid, T can be replaced by 

Fgp[J] = (/ P + /p 2 )/p 3 /p4 - fpfpitfps + /pJ- In the 
context of wave turbulence, the collision integral with 
Tqp is derived based on a statistical averaging, resulting 
in a random phase approximation 1|| 2^]. The collision 
integral in this form was also re-derived recently in a 
microscopic study of the expansion of a matter wave- 
packet based on the GPE [l4| . Another limit of interest 
is the limit of low occupancy (LO), / p C 1, / Pi < 1. In 
this case, F may be replaced by T LO [f] = / P3 /p 4 -/ P / P2 ■ 
Diffusive limit: We are interested in the limit in 
which the average kinetic energy, e, is already the largest 
energy scale. In particular, we will assume that e is much 
larger than the potential energy ~ gn(r). As we will 
see below, under this condition the scattering rate due 
to the inter-particle interaction, l/r co ;, is much smaller 
than £. (In other words, we can speak about states with 
a well defined energy.) Besides, we will assume that 
e is much larger than the disorder induced scattering 
rate l/r p , where p is within the range of typical mo- 
menta which determine e. This inequality not only al- 
lows to develop a comprehensive theory of diffusion but, 
among other things, insures that in dimension d = 2 
the gas cloud expands for a long time before localiza- 
tion effects become essential. Further, we focus on the 
regime for which disorder is the main scattering mech- 
anism, l/r p 3> 1/t co /. It can be expected that under 
these conditions the momentum is randomized on time 
scales of the order of the transport scattering time in- 
duced by the impurities and, therefore, the zeroth an- 
gular harmonic /□ = (f p ) is the dominant one. This 
motivates the standard ansatz used in the transport the- 
ory, / p (r,i) = f (p,r,t) + d nfi(p,r,i) where d is the 
dimensionality, p = |p|, n = p/p and |fi| <C /o- Insert- 
ing this ansatz, a set of two equations can be derived by 
projecting the kinetic equation onto the zeroth and first 
angular harmonic, 

d t fo + «Vfi - V k U (dpf* + ) - (C»[f\)n . 

d t fi + Jv'/o - v*ir - d d P h = ~ + < n ^2 2 [/]) n . (6) 

Neglecting the small dt£[ and (n l C22[/]) n (recall that 
T p is the shortest time-scale), one may resolve for f{ = 
~ ^ l Ud Ep ) fo. Here, and from now on, we as- 
sume that the energy spectrum is quadratic, £ p = p 2 /2m. 
Then, entering the first equation, one gets 

dtfo - (V - VUd £p ) V ep (V - VUd ep + T £p VU) f 
= UM, (7) 
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where D £ = 2eT/(dm) is the diffusion coefficient, T £ — 
^2j^ = — d e hxi>(s), v(e) is the density of states and 
Zdlfo] — C-2 2 [fo}- Indeed, after entering with the ansatz 
for f p into the averaged collision integral (C 2 2[/]) n , the 
largest contribution comes from those terms in which all 
factors of f p are replaced by fo. 

Collision integral: Since fa depends only on the mod- 
ulus of p, but not on n, the angular integrations over 
112, ri3 and ru in the expression for (C22\h\)„ may be 
performed using a method devised in Ref. [21|. T he ex- 
pression for Td [fo] can be written in the form [22j j 

Id[fo] = C d j ' de 2 de^d£ A 8{e p + e 2 - £3- £4) 

x Xd(p,P2,P3,P4)J r [fo}- (8) 

The additional constraint imposed by the energy con- 
servation in the collision integral leads to further simpli- 
fications in the expression for the kernel X3 : 

X3(.P,P2,P3,P4) = y / min(ep,e 2 ,£3,£4)/£p, (9) 



where p t = p £i = ^2/mel. The normalization constant for 
d = 3 is C3 = g 2 m 3 /2ir 3 . In two dimensions, the kernel 
X2 can be found as follows (0 is the Heaviside function): 

X2(P,P2,P3,PA) = P 3+P4-|P-P2| P+P2-|P3-P4| ( 10 ) 

K (A/A) /A, A > A 
A' (A/ A) /A, A < A ' 

where A 2 = PP2P3PA and 

16A 2 = nti (pi +P2 +P3 +P4 - 2jh)\ p1=p . (11) 

The function K is the complete elliptic integral of the first 
kind. The kinematic constraints ensure that each of the 
factors on the RHS in the definition of A 2 is nonnegative. 
The normalization constant is C2 — g\rr? jit . 

Relation to previous results: Let us establish a con- 
nection with the diffusive kinetic equation that was pre- 
viously derived microscopically from the GPE, first in 
the collisionless limit [HI, [H| , and then in the presence 
of collisions [3, 0] • Upon setting n' e = fo (p £ = \j2me) 
and 2TTi/(e)n'(e) = n £ in Eq. ([7]), one immediately arrives 
to the equation 



d t n £ - (V - VUd £ 
= 2nv{e)l d [n'l 



T £ VU)V £ (V - VUd £ )n £ 



(12) 



where the collision integral contains J- as defined in 
Eq. which is valid for arbitrary occupation numbers. 
In the case of large occupation numbers, n'(ei) 3> 1, one 
may substitute T — > Top. The collision integral in this 
form, 2 GP , was derived microscopically in Ref. Q 

As a second example, we make contact with a work 
about thermalization of a rarified gas during a stationary 
diffusion process, see Ref. [lj|. To this end, we specialize 



on I LO , which contains Flo as appropriate for a rarified 
gas with low occupation numbers. For a steady state 
process, in the d = 3 case, the kinetic equation ([7]) can 
be written as 

- (V - VUd e - T £ WU) V £ (V - VUd £ ) n £ 
= 2tw(e) (Z3 [ri; in] - I 3 LO [ri; out]) . (13) 

Here, we split the collision integral explicitly into "in" 
and "out" scattering terms, which take the form 



l£ u [n;m} = -± \ de^de^ /(e; e 3 , e 4 ) n £3 n £i 



n ; out] 



4tt 



de 2 g(e,e 2 ) n £2 n £ 



(14) 



where 



g(e,e 2 ) = / de 3 sj min(e, e 2 , £3, £ + £2 - £3)/£ 2 £2, (15) 



/(e; £3, £4) = V min (e, e 3 + e 4 - e, £3, £4) /££3£4- (16) 

The above formulas are written under the condition that 
all energy arguments are larger than zero. Up to the 
normalization, the functions / and g coincide with the 



functions / and g obtained in Ref. 15( using a diagram- 
matic approach to describe energy transfer processes. 

Mean squared radius: Although the diffusive kinetic 
equation (jTJl constitutes a considerable simplification 
compared to the Boltzmann equation (jTJ) , it is still diffi- 
cult to solve. It turns out, however, that for the special 
case of a two-dimensional system with a constant trans- 
port scattering time and for U = i9 an analytical result 
for the time dependence of the mean squared radius of 
the evolving gas can be found 



d t (r 2 t ) =2dD £ 



(17) 



where e to t — E to t/N is the total energy per particle (for 
the derivation, see 122]). This relation has already been 
noted in Refs. [l2|, 13| f° r theories based on the GPE. 
Here, we generalize this result to the range of applicabil- 
ity of the Boltzmann equation. 

Collision rate: The collision integral contains two 
terms describing the "in"- and "out" -collision channels. 
To estimate the scattering rate 1/t co ii of a state e, let us 
focus on the "out"- term, which correspond to the last 
term in the expression for F[f], see Eq. (j4j. We will be 
interested in the GP-limit of large occupations. Natu- 
rally, the "out"- term is proportional to n £ ; we will write 
it as n £ /T co u. Let us denote the typical kinetic energy 
per particle at point r as e(r). For a conservative es- 
timate of the scattering rate, let us consider an energy 
£ ~ £~(r). Since one has to integrate a product of two dis- 
tribution functions over energies, this ultimately yields a 
factor n 2 (r,t) in the scattering rate: 



1 

Tcoll 



2 n 2 {v,t) 



9d^ 



£(r) 



(18) 
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We consider the case when gd,n{r, t) <C e(r) that leads to 
1/tcoII < 

Discussion: As is well known, for a Bose gas confined 
in a given volume, the redistribution induced by collisions 
leads to a seeding of a state with macroscopically large 
occupation at the bottom of the energy band 0, HH . In 
contrast, in the case of a disordered medium, particles 
that get into states with energies e p < 1/t p as a result 
of collisions are prone to localization. Let us denote the 
energy for which e v = 1/t p as 1/r*. One may try to 
estimate the rate of generation of the density of localized 
particles by integrating the in-scattering term in the col- 
lision integral over the interval of energies e < 1/t*, i.e., 
dnioc/dt &s dn/dt\ e <ci / T «. The analysis shows, however, 
that in 3d, under the conditions chosen for this work, the 
mechanism of seeding of localized particles as a result of 
interparticle collisions is not very effective. 

The situation in restricted dimensions is more involved. 
The difference from 3d is that in d ^ 2 all states are local- 
ized. In 2d the scale of the localization length, h oc {£), in- 
creases exponentially with energy e when it is larger than 
1/t*. In a quasi one-dimensional geometry, as it is well 
known, [ISESl hoc ~ N c hl, where I the disorder scatter- 
ing length at e ~ e, and N c h ~ A/ Xa± ^> 1 is the number 
of channels in a wire or stripe with a cross-section area 
A; a± = min(a z , A). (The strictly Id case will not be dis- 
cussed; a± = a z corresponds to a stripe.) Let us now con- 
sider a Bose-gas cloud of N particles in d ^ 2 which can 
be characterized by a radius R(t) at a time t. Roughly 
speaking, during the course of expansion, the states with 
the energy e are available for localization at the time t s 
when R(t e ) becomes > h oc {s)- This, however, does not 
necessarily mean that such particles are indeed localized 
at this time. For localization, which implies a well defined 
state, it is necessary that the energy uncertainty should 
be smaller than the level spacing A. Hence, for a parti- 
cle with a typical energy e, to start localizing at time tj-, 
the condition l/r co i(t T ) < A /oc = I / [u d lf oc (e)] should be 
fulfilled; here v& is the effective density of states. A nec- 
essary condition for localization of the main population 
of particles into single-particle states can be formulated 
as follows: 



(R 2 (t)) 



> 



max 



(d=l,2)] 



(19) 



where R\ 



.(2d) = ^\l loc (s)N, and R\ (Id) 



(j^) 2 XlN 2 . Thus, if the number of particles is large 
enough, e.g., in Id when N/N ch > (a x /a s )(l/ A) 1 / 2 , the 
expanding cloud will pass hod's) without being stopped. 
The frequent collisions will make the interference, which 
is responsible for localization, ineffective 28 1. 

After reaching R\ , the expansion of the cloud does 
not halt. There are two mechanisms acting in paral- 
lel at t > t-g. One is the expansion of the cloud due 
to the collision-induced diffusion, with the diffusion co- 
efficient D col ~ l 2 loc {e) / t co i ~ l 2 loc (e)g 2 d n 2 {r,t) /e (HH). 



The other mechanism is releasing potential energy stored 
in the cloud which generates a force — Vi?(r, t) per par- 
ticle. As a result the cloud will expand according to a 
power- law rate (i? 2 (i)) ~ t a , with some a < 1. Let us 
discuss a for various mechanisms and different dimen- 
sions. 

Since in the limit of large occupancy D co i oc n 2 (r,t), 
we have to discuss non-linear diffusion described by a 
porous medium equation (PME)[3^, 34 1: 



dn/dt = £>V(n m Vn), 



(20) 



where T> is some constant. The asymptotic solution 
of the PME is well known [H n\r,t) = [(C - 

kr 2 /t a )/t 1 - a ] x / m > 0, where a = 2 /[dm + 2], and C 
and k are constants. We are now interested in m = 2. 
Then the obtained expression for a coincides with the ex- 
ponent a = l/(d+l) conjectured in Ref. 1 3 71 for spreading 
of a wave-packet in non-linear disordered lattices in the 
case when a mechanism of dephasing has been introduced 
;38]. In d = 2 it yields a = 1/3. 

To get more insight, let us discuss a Bose gas of large 
occupancy diffusing in a broad stripe. Such a system cor- 
responds to the Id GPE/NLSE. However, it is natural to 
expect that the presence of a large number of transverse 
channels will lead to dephasing. For this case, we get 
q = 1/2, which is in full correspondence with the expo- 
nent a observed numerically in Ref. 1371 . Furthermore, we 
can compare this exponent with the rate of spreading of 
a cloud of the Bose gas caused by the force — Vi9. This 
yields a current density j(r,t) oc =D co ;V?9(r, t) which 
corresponds to m = 3 and leads to a = 2/5. This ex- 
ponent coincides precisely with the one conjectured by 
Shepelyansky [39|. We find it natural, because his ar- 
guments where concentrated on the energy balance in a 
spreading packet. 

A technique of deriving the collision integral in clas- 
sical equations (NLSE/GPE) was described in our pre- 
vious paper [l4j |. We believe that an analysis in terms 
of collisional kinetics helps to clarify the physics of sub- 
diffusive spreading of wave packets in disordered non- 
linear lattices. Quite generally, one may guess that the 
collision-induced diffusion is responsible for the spread- 
ing. A properly designed experiment, with a large enough 
number of Bose-atoms and a large number of channels 
may help to get valuable information about transport 
properties of disordered interacting systems, both quan- 
tum and classical. 
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SUPPLEMENTAL MATERIAL 
Derivation of relation ([8]) for the collision integral T c i. 

In this section, we provide some additional informa- 
tion on the derivation of Eq. ([5]) for the collision integral 
^d[/o] = C| 2 [/o]. Starting point is Eq. d3j. After replac- 
ing f —> fo the angular averages in n 2 , 113 and 114 can be 
performed. To this end let us define 

R d = J dn 2 rfn 3 dn 4 S(p + p 2 - P3 - Pa)- (21) 



After rewriting 



s(p + P2 p 3 - p 4 ) = J (d\y x ^-^-^\ 



(22) 



6 



one finds for the angular average 
Rd = fid 



-fTT^-Vd ( Ap)r? d (Xp 2 )Vd(Xp3 )Vd(Xp 4 ) ■ (23) 
o l /7r J 



Here, fid is the solid angle in d dimensions, and r)d(Xp) 
(e ±lAp ) n . Further, 



sm Ix) 



772(a) = J (x), 



(24) 



where Jo is a Bessel function of the first kind [401 ] ■ In 
three dimensions, the remaining integral in A is elemen- 
tary, while for two dimensions it can be found in Ref. |41 
This brings us directly to Eq. ([8]) with the kernels \d 
given by Eqs. © and (fT0|) . 



Derivation of relation (|17[) for the mean squared 
radius. 



In this section, we present the derivation of relation 
(|17p . which describes the time evolution of the mean ra- 
dius squared for the 2rf-Bose gas with energy-independent 
t. The reasoning is based on the conservation laws for 
particle number and energy. It is crucial for the argu- 
ments to consider the self-induced potential only, U = ■&. 
The continuity equation for the density is easily obtained 



from the Boltzmann equation upon integration in p: 



d t n + Vj = 0, 



(25) 



where 



n = (/p) p and J = (/p v ) p — V/p 

J(dp)(. ..). To obtain the continuity equation for the 
energy density, 



and 



. ) means 



dtPE + Vjb = 0, 



(26) 



one multiplies the Boltzmann equation by e p and inte- 
grates in p. It gives 

PE = (e P f P ) p +g d n 2 , 3e - (ve p / p ) p + tfj. (27) 

We immediately conclude that the total energy E to t = 
J dvpE is conserved. The mean radius squared is defined 
as (r|) = jj J drr 2 n(r). From the continuity equation 
for the density it follows that d t {r 2 ) = jj f drrj(r,t). 
Now consider the mean radius squared in the diffusive 
limit, for which j = (vfi) p . Two conditions are needed 
for the argument: (i) T e = 0, i.e., the density of states 
is constant, and (ii) D e = ke, where k is a constant. 
Under these conditions j = — k\7pE- Plugging this into 
the relation for (rf) one immediately finds relation (fP7j) 
with etot = Etot/N. Both conditions (i) and (ii) are 
fulfilled for the 2d-case with energy-independent r. 



